Krylov subspaces and the analytic grade

Typical behaviour of the solution of a linear system of equations obtained iteratively by Krylov methods can be characterized by three stages. Initially the residual diminishes steadily; this is followed by stagnation and finally rapid convergence near the algebraic grade. This study examines this behaviour in terms of the concepts of approximately invariant subspace and what we have called the analytic grade of a Krylov sequence. It is shown how the small Ritz values play a vital role in the convergence and how this knowledge helps in the construction of an effective preconditioner. Copyright © 2004 John Wiley & Sons, Ltd.     

Adaptive remeshing for ductile fracture prediction in metal formingRemaillage adaptatif pour la prévision de la rupture ductile en mise en forme

The analysis of mechanical structures using the Finite Element Method in the framework of large elastoplastic strain, needs frequent remeshing of the deformed domain during computation. Indeed, the remeshing is due to the large geometrical distortion of finite elements and the adaptation to the physical behavior of the solution. This paper gives the necessary steps to remesh a mechanical structure during large elastoplastic deformations with damage. An important part of this process is constituted by geometrical and physical error estimates. The proposed method is integrated in a computational environment using the ABAQUS/Explicit solver and the BL2D-V2 adaptive mesher. To cite this article: H. Borouchaki et al., C. R. Mecanique 330 (2002) 709–716.     

用模拟万用表测电容

介绍了用模拟万用表和秒表测量电容的简单方法，对测量结果的误差原因进行了分析，给出了提高测量精度的几点措施。     

Assessing the Error Probability of the Model Selection Test

The asymptotic error probability of Linhart's model selection test isevaluated, and compared with the nominal significance level. We examine thecase where the expected discrepancies of the candidate models from the truemodel are asymptotically equal. The local alternatives method is employed inthe limiting operation of the asymptotic evaluation. Although the errorprobability under the null hypothesis is actually shown to be equal to orless than the level for most situations, intolerable violations of the errorcontrol are observed for nested models: It is often erroneously concludedthat the smaller model is significantly better than the larger model. Toprevent this violation, a modification of Linhart's test statistic isproposed. The effectiveness of the proposed test is confirmed throughtheoretical analysis and numerical simulations.     

一类常微分方程组的紧支撑样条小波方法探讨

应用三次紧支撑样条小波插值函数得到了求一类常微分方程组数值解的隐式公式,并求得到其局部截断误差为O(-h5).在此基础上给出1个显式校正求解公式,并讨论得到其局部截断误差为O(-h4).     

A mixed‐FEM formulation for nonlinear incompressible elasticity in the plane

This article deals with an expanded mixed finite element formulation, based on the Hu‐Washizu principle, for a nonlinear incompressible material in the plane. We follow our related previous works and introduce both the stress and the strain tensors as further unknowns, which yields a two‐fold saddle point operator equation as the corresponding variational formulation. A slight generalization of the classical Babu?ka‐Brezzi's theory is applied to prove unique solvability of the continuous and discrete formulations, and to derive the corresponding a priori error analysis. An extension of the well‐known PEERS space is used to define an stable associated Galerkin scheme. Finally, we provide an a posteriori error analysis based on the classical Bank‐Weiser approach. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 105–128, 2002     

EFFICIENT SIXTH ORDER P-STABLE METHODS WITH MINIMAL LOCAL TRUNCATION ERROR FOR y"=f(x,y)

1. IntroductionWe consider a class of direct hybrid methods proposed in [11 for solving the second orderinitial value problemy" = f(t,y), y(0),y'(0) given (1.1)The basic method has the formandHere t. = nh and we define t.l.. = t. I aih, i = 1, 2 and n=0,1…     

Theoretical and numerical analysis on multispectral bioluminescence tomography

Recently, molecular imaging has been rapidly developed to studyphysiological and pathological processes in vivo at the cellularand molecular levels. Among molecular imaging modalities, opticalimaging has attracted a major attention for its unique advantages.In this paper, we establish a mathematical framework for multispectralbioluminescence tomography (BLT) that allows simultaneous studiesof multiple optical reporters. We show solution existence, uniquenessand continuous dependence on data as well as the limiting behaviourswhen the regularization parameter approaches zero or when thepenalty parameter approaches infinity. Then, we propose twonumerical schemes for multispectral BLT and derive error estimatesfor the corresponding solutions.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

A lower bound on the performance of the quadratic discriminant function

The quadratic discriminant function is often used to separate two classes of points in a multidimensional space. When the two classes are normally distributed, this results in the optimum separation. In some cases however, the assumption of normality is a poor one and the classification error is increased. The current paper derives an upper bound for the classification error due to a quadratic decision surface. The bound is strict when the class means and covariances and the quadratic discriminant surface satisfy certain specified symmetry conditions.